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ABSTRACT 


Quick-Mixing Studies Under Reacting Conditions 

The low-NO x emitting potential of rich-burn/ quick-mix/lean-burn 
(RQL) combustion makes it an attractive option for engines of future 
stratospheric aircraft. Because NO x formation is exponentially dependent on 
temperature, the success of the RQL combustor depends on minimizing high 
temperature stoichiometric pocket formation in the quick-mixing section. An 
experiment was designed and built, and tests were performed to characterize 
reaction and mixing properties of jets issuing from round orifices into a hot, 
fuel-rich crossflow confined in a cylindrical duct. The reactor operates on 
propane and presents a uniform, non-swirling mixture to the mixing modules. 
Modules consisting of round orifice configurations of 8, 9, 10, 12, 14, and 18 
holes were evaluated at a momentum-flux ratio of 57 and jet-to-mainstream 
mass-flow ratio of 2.5. Temperatures and concentrations of 0 2 , C0 2 , CO, HC, 
and NO x were obtained upstream, downstream, and within the orifice plane to 
determine jet penetration as well as reaction processes. Jet penetration was a 
function of the number of orifices and affected the mixing in the reacting 
system. Of the six configurations tested, the 14-hole module produced jet 
penetration close to the module half-radius and yielded the best mixing and 
most complete combustion at a plane one duct diameter from the orifice 
leading edge. The results reveal that substantial reaction and heat release 


occur in the jet mixing zone when the entering effluent is hot and rich, and 
that the experiment as designed will serve to explore satisfactorily jet mixing 
behavior under realistic reacting conditions in future studies. 
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C d discharge coefficient 

DR jet-to-mainstream density ratio 

FAR fuel-air ratio 

d orifice axial height, or round hole diameter 

J jet-to-mainstream momentum-flux ratio = (pV 2 )- et / (pV 2 ) main 

MR jet-to-mainstream mass-flow ratio 

P pressure 

ppm part per million by volume 

r radial distance from center of module 

R radius of the quick-mix module 

R 1 radius one-third of R from the center 

R, radius two-thirds of R from the center 

SLPM standard liters per minute 

Tj temperature of component i (i = jet, mainflow) 

V ref reference velocity 

VR jet-to-mainstream velocity ratio 

z axial distance from leading edge of orifice 

<f) equivalence ratio = ( fuel / air) actua , / (fuel / air) stoichiometnc 
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CHAPTER 1 


INTRODUCTION 


1.1 Overview 

Jet mixing with a crossflow has been well-researched because of its 
occurrence in a wide range of applications such as pollution control, vertical 
short takeoff/landing (VSTOL) aircraft, and gas turbine design. 
Understanding the dynamics of jets interacting with a crossflow is essential in 
predicting jet trajectory and mixing in these systems. 

The gas turbine combustor relies heavily on jet-crossflow mixing to 
achieve flame stability, completion of reaction, and cool uniform exit 
temperatures. The importance of jet mixing with a crossflow is further evident 
in the Rich burn-Quick mix-Lean burn (RQL) combustion concept. The RQL 
scheme is a low-NO x combustor being considered for powering the next 
generation fleet of supersonic aircraft (Shaw, 1991). The premise behind the 
RQL combustor lies in staging the combustion process in fuel-rich and fuel-lean 
zones to avoid the high temperatures associated with near-stoichiometric 
combustion. High temperatures encourage the production of NO x , a pollutant 
which when released into the upper atmosphere participates in the destruction 
of stratospheric ozone. 

The success of the RQL combustor rests with the performance of the 
quick-mixing section that bridges the rich and lean zones. Rapid and thorough 
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mixing of jet air with a crossflow of rich products to complete the reaction is 
desired in order to decrease the formation and residence time of stoichiometric 
pockets of fluid. These fluid pockets are undesirable because NO x production 
is accelerated by the high temperatures induced by stoichiometric conditions. 

The key to reducing NO x formation in the quick-mixing section lies in 
determining the effect of orifice configurations on jet penetration and mixing 
uniformity. Investigations have focused on varying flow and geometric 
parameters to determine configurations that lead to uniform mixing within a 
specified duct length. Multiple jet mixing experiments have mainly been 
performed under non-reacting rather than reacting conditions because fewer 
complications are involved. Reacting flow investigations reported in the 
literature have mainly consisted of numerical studies, and the few existing 
experimental reacting studies have not been directed at understanding the 
dynamics of jet mixing with a rich reacting crossflow. 

1.2 Research Goal and Objectives 

Previous studies conducted by Hatch, et al. (1996), Kroll, et al. (1996), and 
Sowa, et al. (1994) investigated the mechanistic properties of jet mixing with a 
heated cylindrical crossflow under non-reacting conditions. The current 
research project builds upon the non-reacting experiments by initiating 
performance studies of the cylindrical quick-mixing region under reacting 
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conditions, with the goal being to characterize jet mixing and reaction for 
selected orifice configurations. 

The objectives that need to be accomplished in order to attain the 
research goal are: 

1) Conduct a detailed literature review on research and issues related to the 

RQL combustor. 

2) Design, construct, and validate a test stand with reacting flow capability. 

3) Ensure experimental apparatus integrity and rich product uniformity. 

4) Design a test protocol and matrix, and conduct reacting flow experiments to 

measure temperature and species concentration profiles. 

5) Analyze acquired data. 


* 


CHAPTER 2 


BACKGROUND 

2.1 High Speed Civil Transport Program 

A joint effort between the National Aeronautics and Space 
Administration (NASA) and aircraft-related industries has been formed to 
develop technology needed to support production of a fleet of High Speed 
Civil Transport (HSCT) aircraft by the twenty-first century. This second- 
generation civilian supersonic aircraft (Figure 2.1) is planned to become more 
commercially successful than its Concorde predecessor with its improved fuel 
efficiency and capability to transport three times as many passengers at twice 
the distance (Strack and Morris, 1988). Certain technological issues need to be 
resolved, however, before such a fleet of supersonic aircraft can be realized. 



Figure 2.1 High Speed Civil Transport 
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The High Speed Research Program (HSRP) was initiated in 1990 by 
NASA to study the issues involved with developing and supporting a fleet of 
advanced supersonic aircraft. Although technical and economic issues require 
consideration in ensuring the viability of the program, environmental barriers 
relating to noise and emissions are preventing the aircraft from leaving the 
ground. While the noise generated from the HSCT aircraft mainly impacts 
communities in the flight path or vicinity of airports, emissions into the 
atmosphere may produce adverse effects over a global area and population. 

The HSCT aircraft is designed to cruise between 18.3-27.4 km (60-90,000 
ft.) at speeds of Mach 2-3 for optimal fuel economy (Koff, 1994). The planned 
flight altitude falls in the domain of the stratosphere, a stable region in the 
atmosphere where the ozone layer resides. In the 1970s an evaluation of the 
impact of stratospheric flight on the environment was carried out by the 
Climatic Impact Committee formed by the National Research Council, the 
National Academy of Sciences, and the National Academy of Engineering 
(Grobecker, et ah, 1974). The committee determined that of all the emissions 
from stratospheric aircraft (including water, carbon dioxide, carbon monoxide, 
sulfur oxides, nitrogen oxides, and soot), nitrogen oxides (NO x ) were a primary 
concern because of their deleterious effect on the stratospheric ozone layer 
(Climatic Impact Committee, 1975). Due to current adverse public opinion 
toward environmentally detrimental technologies, the industrial sector is 
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hesitating in developing the HSCT aircraft until its potential for depleting 
stratospheric ozone is resolved. 

2.2 Stratospheric Ozone Depletion 

The stratosphere contains up to 90% of the ozone (0 3 ) in the entire 
atmosphere, with maximum concentrations found between 15-25 km in the 
lower stratosphere (Rowland, 1991). The ozone that resides in the stratosphere 
forms a thin blanket over the Earth which absorbs the UV-B waveband, or 
ultraviolet radiation shorter than 320 nm. UV-B radiation destroys cells of 
plants and animals and leads to skin cancer, eye cataracts, and deterioration of 
the immune system in humans. The destruction of the ozone layer increases 
terrestrial exposure to harmful UV-B radiation and increases urban air 
pollution caused by the photolysis of formaldehyde in photochemical smog 
(Masters, 1991). On an atmospheric scale, the gradual loss of the ozone layer 
decreases the stratospheric temperature and consequently, the circulation in the 
atmosphere (Wayne, 1985). 

The stratospheric ozone production and destruction cycle can be 
described in the following series of steps (equations 2.1-2.4) proposed by 
Chapman (1930): 


0 2 + UV (k < 242nm) O + O 

(2.1) 

0 + 0 2 + M-»0 3 + M 

(2.2) 

o 3 + UV -> 0 + o 2 

(2.3) 
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O + Oj— > 0 2 + 0 2 (2.4) 

where X represents the radiation wavelength. The Chapman reactions, 
however, consistently overpredict the amount of ozone recorded in 
experiments. The explanation for the occurrence of lower than expected 
concentrations of ozone was first suggested by Bates and Nicolet (1950) as 
being caused by destructive catalytic processes. The following series of free 
radical catalytic reactions (equations 2.5-2. 6) reflecting the additional 
destruction of ozone was then appended to the sequence to compensate for the 
overprediction: 


x + o 3 -» xo + o 2 

(2.5) 

xo + 0 x + o 2 

(2.6) 

O + 0 3 — > 0 2 + 0 2 

(2.7) 


where X is a catalyst representing either HO x , NO x , or C10 x . The pair of ozone- 
destroying reactions in equations 2.5 and 2.6 can continue indefinitely until the 
catalyst X is removed. Unfortunately, the catalyst may linger for a long period 
of time in the stratosphere due to the temperature inversion which inhibits the 
vertical mixing that would normally aid in its removal. 

The role of NO x in stratospheric ozone destruction via the catalytic 
reactions emphasizes the need to reduce NO x emissions from the gas turbine 
combustors that will power the HSCT aircraft. In order to control NO x 
production during the combustion process, the mechanism and chemistry 
behind NO x formation must be understood. 
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2.3 NO x Formation 

Oxides of nitrogen (NO x ) are pollutants formed during the combustion 
process that are partially responsible for the degradation of the ozone layer. 
The NO x label represents both nitric oxide (NO) and nitrogen dioxide (N0 2 ). 
However, because N0 2 comprises less than 10% of the total amount of NO x , 
NO x is almost all NO. 

During the combustion process, NO x can be formed via three pathways 
known as prompt, fuel, and thermal NO x . Prompt NO x is formed in the 
primary reaction zone during the initial stages of combustion when 
hydrocarbon radicals or fuel fragments attack atmospheric nitrogen molecules 
(N 2 ) (Fenimore, 1971). The resulting atomic nitrogen product from the 
dissociation of N 2 reacts with oxygen molecules to form NO. Prompt NO x 
formation, however, comprises a relatively small fraction of the total NO x 
formed. 

Fuel NO x is formed when the nitrogen in the fuel is oxidized. The fuel- 
bound nitrogen is typically bonded to carbon and hydrogen in the form of 
ammonia, pyridine, and other amines (Glassman, 1987). Fuel NO x can be 
controlled by choosing fuels with lower nitrogen content. 

Thermal NO x formed from atmospheric nitrogen is the main source of 
NO x emissions for combustion systems operating at high temperatures and 
with long residence times. The series of reactions listed as equations 2.8 and 
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2.9, known as the Zeldovich mechanism (Zeldovich, 1946), describes the 
production of thermal NO x . 


N 2 + O NO + N 

(2.8) 

N + 0 2 NO + O 

(2.9) 


The first reaction (equation 2.8) is the rate-limiting step because it requires a 
high activation energy to initiate the reaction by breaking the triple bond that 
holds the nitrogen molecule together. 

For rich and near-stoichiometric flames the second reaction (equation 
2.9) can be replaced by the extended Zeldovich mechanism listed as equation 
2 . 10 . 

N + OH — » NO + H (2.10) 

The dependence of the NO reaction rate on temperature is shown in 
Figure 2.2. Above a temperature of 1800K (2780°F) the reaction rate constant 
and hence production of NO increases exponentially. To control thermal NO 
production, the temperature of the reaction should be maintained below this 
critical temperature. 

The reaction temperature is dependent on the fuel-air equivalence ratio 
<J), defined in equation 2.11 as 

<t> = (fuel/air),^ (2.11) 

(fuel / a i r ) s t 0 i C i 1 i ornetr ic 
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Temperature (K) 

Figure 2.2 Effect of Temperature on NO Reaction Rate 
(Adapted from Samuelsen, 1975) 

or the overall ratio of the actual to stoichiometric fuel-air ratios. A reaction that 
is stoichiometric (<t> = 1) contains no excess fuel or air in the combustion 
products. At near-stoichiometric conditions the adiabatic flame temperature, 
or the highest theoretical reaction temperature, is attained. As the reactant 
composition approaches fuel-rich (((> > 1) or fuel-lean (<f> < 1) conditions, 
temperatures decrease sharply. The NO formation dependency on equivalence 
ratio, shown in Figure 2.3, shows a bell-shaped curve that also corresponds to 
flame temperature dependency on equivalence ratio. The peak of the curve is 
shifted slightly toward the lean condition due to an abundance of oxygen 
radicals. This O-atom overshoot slightly increases the NO formation rate at 


leaner conditions. 
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Figure 2.3 Effect of Equivalence Ratio on NO Formation 
(Adapted from Samuelsen, 1975) 

Because high inlet temperatures from the compressor stage contribute to 
the even higher combustion temperatures occurring in the combustor, the 
thermal mechanism is the primary mode of NO x production in a gas turbine. 
Control of reaction temperature by varying the equivalence ratio and 
minimizing the residence time of the reaction at high temperatures is the key to 
reducing NO x production rates in a gas turbine combustor. 

2.4 The Gas Turbine Combustor 
2.4.1 Description 

In a typical aircraft engine the combustor is situated between the 
compressor and turbine (Figure 2.4). Using the hot air from the compressor 
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stage, the combustor transforms the chemical energy in fuel into the heat 
energy that drives the turbine. The basic gas turbine combustor (Figure 2.5) 
consists of primary, intermediate, and dilution zones. The primary zone 
houses the dome region in which liquid fuel is vaporized and mixed with air. 
Additional air is added to the partially combusted products in the intermediate 
zone to complete the reaction. The dilution zone mixes in air with the complete 
combustion products to tailor the gas temperatures to turbine blade material 
specifications. 



Figure 2.4 Aircraft Engine Schematic (Adapted from Koff, 1994) 
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Figure 2.5 Combustor Parts and Zones (Adapted from Lefebvre, 1983) 


The combustor operates via a process of continuous combustion in 
which fresh reactants are continually injected, burned, and exhausted from the 
combustor (Figure 2.6). A flame stabilizer such as a bluff body or a swirler 
creates a recirculation region in the primary zone that funnels hot reactants 
back to the dome region. The hot reactants provide a source of ignition for the 
fresh fuel and air mixture injected into the combustor. 
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Air Air 



Figure 2.6 Continuous Combustion Aerodynamics 
(Adapted from Samuelsen, 1975) 


2.4.2 Ultra Low-NO x Combustors 

Current commercial engines produce a NO x Emissions Index (El), which 
is a mass-based measurement of pollutant emitted for a given amount of fuel, 
between 40-60 g NO x /kg fuel. Based on this El, a projected fleet of 500 
supersonic aircraft flying in the lower stratosphere could decrease ozone levels 
annually by 20% on a global basis. However, if the HSCT aircraft complies 
with the HSRP goal of an El of 5 g NO x /kg fuel, the potential destruction of 
ozone could be as low as 2-3% (Johnston, et al, 1989). 
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To attain the HSRP El goal without sacrificing engine efficiency, ultra 
low-NO x combustor technologies are being investigated. The most promising 
concepts for meeting the El requirement are the Lean Premixed-Prevaporized 
and the Rich burn-Quick mix-Lean bum systems. These ultra low-NO x 
concepts reduce thermal NO x formation by operating at lean or rich 
environments to take advantage of the lower temperatures associated with non- 
stoichiometric equivalence ratios. 

The Lean Premixed-Prevaporized (LPP) concept (Figure 2.7) involves a 
single stage of combustion under fuel-lean conditions. The liquid fuel is 
completely vaporized and mixed thoroughly with air before combustion. NO x 
formation is reduced since the lean operating condition reduces flame 
temperatures while the avoidance of droplet burning eliminates near- 
stoichiometric fluid pockets. The LPP concept, however, faces technological 
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Figure 2.7 Lean Premixed-Prevaporized Concept (Adapted from Shaw, 1991) 


16 


and safety issues such as incomplete fuel-air mixing, the potential for 
autoignition due to high inlet temperatures, flashback (flame propagation back 
to the fuel source), flame blowout at low power conditions, and poor re-light 
capability (Lefebvre, 1983). 

The Rich burn-Quick mix-Lean burn (RQL) combustor (Figure 2.8) has a 
wider combustion stability limit and does not incur the hazards of autoignition 
apparent in the LPP combustor. The RQL combustor was originally conceived 
as a means to control fuel NO x because a fuel-rich environment decreases fuel- 
bound nitrogen conversion to NO x (Tacina, 1990). However, the concept is 
being applied to thermal NO x reduction because it operates on a two-stage 
combustion process which limits the time spent at near-stoichiometric 
conditions. Fuel and air are first burned in a fuel-rich environment with a 
limited amount of oxygen available for NO production. The rich products are 
then rapidly mixed with jets of air in the quick-mixing region to bring the 


Quick-Mixing Jets 



Rich Product 

► 


O 

O 


Lean Product 
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Figure 2.8 Rich Burn-Quick Mix-Lean Burn Concept 
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reaction to completion in a fuel-lean environment. The quick-mixing section is 
necked-down to prevent the backflow of products into the rich zone. 

In a literature review of experiments performed on the LPP and RQL 
combustors, Tacina (1990) determined that the LPP combustor produces a 
lower amount of NO x than the RQL combustor. The higher NO x emissions 
from the RQL combustor was attributed to the near-stoichiometric fluid pockets 
formed in the quick-mixing section. Though the RQL concept still possesses an 
ultra low-NO x emission potential, its ultimate success hinges on the 
performance of the quick-mixing section that bridges the rich and lean zones. 
The mixing of jet air with a rich effluent must be performed rapidly and 
uniformly in order to decrease the time of transition between the two zones 
and the occurrence of near-stoichiometric fluid pockets. The challenge of 
optimizing the mixing process in the quick-mix section of the combustor rests 
with understanding the mechanism behind jet mixing in a crossflow. 

2.5 Jets in Crossflow 

2.5.1 Single Jet Structure 

The jet-in-crossflow problem has been studied extensively due to its 
broad range of applicability to such diverse fields as gas turbine cooling and 
staging, fuel-air premixing, vertical short takeoff/ landing (VSTOL) aircraft, 
and pollutant discharge from stacks or pipes. 
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A round jet entering a crossflow forms a complex three-dimensional free 
turbulent shear flow (Figure 2.9). The structure of the jet in a crossflow is 
dependent on the interplay between the jet and crossflow momenta. Toward 
the base of the jet, the crossflow splits around the jet flow as it would around a 
solid body. Farther up along the jet, mixing occurs between the two fluids. 


Velocity 



Figure 2.9 Round Jet Deflection and Crossflow Entrainment 
(Adapted from Lefebvre, 1983) 


A reverse flow region occurring in the wake accelerates fluid from the 
crossflow and entrains the fluid into the jet. The vorticity of the crossflow 
interacts with the circular vortical field of the jet to produce a bound vortex 
shaped like a horseshoe. The vortical field induces the formation of a pair of 
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counter-rotating vortices within the jet that deform the round jet into a kidney 
shape (Andreopoulos and Rodi, 1984). Downstream from the jet entrance, the 
diffusion of vorticity weakens the pair of vortices: the jet loses its coherency 
and subsequently disperses into the crossflow. 

The characterization of the single jet in a non-reacting crossflow has 
been pursued in experimental and analytical modeling studies. Turbulence 
measurements verifying the structure of the jet-crossflow interaction 
(Andreopoulos and Rodi, 1984; Andreopoulos, 1985; Fearn and Weston, 1974; 
Kamotani and Greber, 1972; Sherif and Pletcher, 1989) as well as temperature 
distributions of heated jet mixing with a crossflow have been obtained 
(Kamotani and Greber, 1972). Analytical models have been developed that 
characterize the vortex behavior of a hot or cold jet in a crossflow (Karagozian, 
et ah, 1986) and that predict the jet vortex trajectory (Karagozian, 1986). These 
theoretical studies contribute to the understanding of the general structure of a 
single jet in a crossflow, and subsequently help in the analysis of multiple jet 
mixing in a crossflow. 

2.5.2 Multiple Jets in a Confined Non-Reacting Crossflow 

The quick-mixing section of the RQL combustor is comprised of a row of 
jet orifices spaced around the circumference of the combustor. The crossflow is 
confined to either a cylindrical or annular combustor geometry. In a confined 
crossflow problem, flow properties such as the jet-to-mainstream density. 
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mass-flow, and momentum-flux ratios as well as the geometries constraining 
the jet and crossflow influence the degree of mixing that occurs. The most 
important flow variable in the confined subsonic crossflow system is the jet-to- 
mainstream momentum-flux ratio J (Holdeman, 1993), defined in equation 2.12 


as 


J 


(pvu 


( 2 . 12 ). 


mam 


The momentum-flux ratio must be determined before an orifice configuration 
of a certain number, shape, and placement can be designed. 

Extensive experimental and numerical studies on jets in a confined 
crossflow have been performed under non-reacting conditions to examine the 
effect of jet orifice configurations on mixing in different duct geometries and at 
various momentum-flux ratios. Non-reacting studies of jets issuing into 
rectangular (Bain, et ah, 1992, 1994, 1995; Liscinsky, et ah, 1992, 1993, 1994, 1996) 
and cylindrical duct geometries (Hatch, et ah, 1995a, 1995b; Howe, et al., 1991; 
Kroll, et ah, 1993; Oechsle, et ah, 1992, 1993; Smith, et ah, 1991; Sowa, et ah, 1994; 
Talpallikar, et ah, 1992; Vranos, et ah, 1991) have been studied for their 
applicability to annular and can combustor configurations. 

Among the primary goals of non-reacting research on jet mixing in a 
confined crossflow is to determine orifice configurations that lead to optimal 
mixing within a specified duct length. In the cylindrical duct geometry, 
experimental surveys of the effect of momentum-flux ratio and the shape. 
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orientation, and number of orifices on mixing were performed by Hatch, et al. 
(1995a) and Kroll, et al. (1993) in order to gain a mechanistic understanding of 
jet penetration and mixing dynamics. A systematic optimization scheme on 

#• 

experimental data was then undertaken by Sowa, et al. (1994) to determine the 

# 

orifice configurations leading to optimal mixing at a set momentum flux ratio. 
A non-linear relationship between orifice shape, number, and orientation was 
revealed with respect to mixing, and allowed for the possibility of more than 
one optimal orifice combination. 

Non-reacting flow experiments have been conducted in lieu of reacting 
experiments in order to benefit from the advantages (less complicated, more 
amenable to diagnostic interrogation, more amenable to modeling) of the non- 
reacting environment. A numerical study by Oechsle, et al. (1994) showed that 
qualitatively similar mixture non-uniformity flow fields were obtained in 
reacting and non-reacting simulations. Another numerical study by 
Talpallikar, et al. (1992) showed non-reacting and reacting flows exhibiting 
optimum mixing at the same momentum-flux ratio for a particular slotted 
orifice configuration. These studies lend credence to the use of non-reacting 
tests as a screening tool for potential RQL mixing configurations. The 
screening potential of the non-reacting tests has led to the development of 
numerical codes that predict NO x production at actual flight conditions based 
on non-reacting mixing parameters (Hatch, et al . , 1995b). Despite the insight 
gained from non-reacting tests on jet mixing, an experimental correlation 
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between non-reacting and reacting tests has yet to be established conclusively. 
Extensive experimental data from reacting tests are required to help validate 
the use of non-reacting tests in predicting mixing under reacting conditions. 

2.5.3 Multiple Jets in a Confined Reacting Crossflow 

Numerical studies on jet mixing in a reacting crossflow have been 
undertaken to characterize the flowfield and NO x production and to relate the 
results to non-reacting flows. Howe, et al. (1991) varied the jet momentum-flux 
ratio in a study comparing non-reacting and reacting cases. In both 
environments, the momentum-flux ratio affected jet penetration depth. The 
reacting case produced jets with a lower penetration depth than the non- 
reacting case because the increase in mainstream velocity from the reaction heat 
release decreased the momentum-flux ratio. 

Oechsle, et al. (1994) found that at set momentum-flux, mass-flow, and 
density ratios, the reacting flow exhibited a lower degree of mixing than the 
non-reacting case. Jet core diffusion and mixing with the crossflow was also 
found not to be as great in the reacting case as in the non-reacting case. The 
studies by Howe, et al. and Oechsle, et al. verify the importance of the 
momentum-flux ratio mixing parameter in affecting jet penetration and the 
degree of mixing. 

Oechsle and Holdeman (1995) performed a numerical reacting flow 
study at HSCT flight conditions in a cylindrical geometry. The momentum- 
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flux ratio and orifice shape were varied, and non-uniformity mixing 
parameters based on equivalence ratio and statistical analysis were used to 
evaluate jet mixing. Results showed that jet penetration affected both the 
mixing flow field and NO x production, as the over- and underpenetration of 
jets led to higher NO x production. 

A few experimental studies have been performed to characterize the 
reacting flow in a model gas turbine combustor (Noyce, et al, 1981; Heitor and 
Whitelaw, 1986). The results, though, are not applicable to the quick-mixing 
regime of the RQL combustor where rapid jet mixing and high momentum-flux 
ratios occur. As this moment, only one experimental study specific to the RQL 
combustor is reported in the open literature. Zarzalis, et al. (1992) performed a 
reacting experiment to determine the effect of different inlet pressures and 
temperatures on NO x emissions. However, their study did not address the 
mechanistic processes governing jet mixing in the quick-mixing section. 

This study builds upon the research of Hatch, et al. (1996), Kroll, et al. 
(1996), and Sowa, et al. (1994) on non-reacting jet mixing in a cylindrical 
crossflow by initiating similar performance studies under reacting conditions. 
The purpose of the study was to obtain, for a rich reacting flow in a cylindrical 
RQL simulation, species concentration and temperature distributions in order 
to evaluate jet penetration and mixing, and provide an initial database for 


numerical simulations. 


CHAPTER 3 


APPROACH 


The goal of characterizing jet mixing in a rich reacting crossflow was 
addressed in four phases: (1) reacting flow facility construction, (2) rich 

product uniformity evaluation, (3) test matrix specification, and (4) data 
analysis. 

Phase 1: Facility Construction 

This phase encompassed the retrofitting of the non-reacting facility used 
previously by Hatch, et al (1996) and Kroll, et al. (1996) for their mixing studies. 
The upgraded facility features a refractory-lined can combustor with a 
removable top section to insert a flow conditioner, an aluminum cylindrical 
chamber that serves as a plenum feed for the quick-mix jets, and quartz 
modules to house the reaction in the quick-mixing regime. 

Phase 2 : Rich Product Uniformity Evaluation 

The design of the rich product generator was an iterative process that 
involved reconfiguring and interchanging the system with various parts to 
produce an experiment that was safe to operate as well as able to produce a 
uniform rich combustion product. 
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Several options were pursued in achieving an ultimate design that 
produced a rich zone equivalence ratio <j) approaching 1.7. Space limitations 

led to the initial use of a Lean Burn Injector (LBI) to mix fuel and air rapidly in 

t 

a short combustor length. However, because the swirling inflow of rich 

S 

products induced by the LBI/swirler assembly was not desirable for the 
objectives of this study, an alternate means of generating rich products was 
pursued. 

A premixed fuel and air system that used a ceramic foam matrix as a 
flameholder for the reaction was proposed to replace the swirler as a flame 
stabilizer. However, the ceramic foam failed to hold the reaction in its porous 
structure at the desired rich equivalence ratio. The resulting product from the 
system was also green, which signified the presence of C 2 radicals in the 
reaction (Glassman, 1987). Species concentration profiles obtained in the green- 
tinged effluent showed high 0 2 concentrations when 0% was expected, high 
unburned hydrocarbon concentrations, and lower than expected CO and C0 2 
concentrations. These observations suggested that the reaction did not attain a 
residence time that was long enough to convert the radicals to the expected rich 
product concentrations. 

One solution to increasing the residence time was to increase the length 

« 

of the combustor, but space limitations precluded this option. The final design 
re-incorporated the swirler in the premixed fuel-air and ceramic foam system. 
Instead of functioning as a flameholder, the foam was used as a swirl 
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dissipater. Several iterations were required to determine the placement of the 
foam with respect to the swirler. The resulting scheme placed the foam five 
duct diameters downstream from the swirler to allow enough distance for the 
recirculation zone. Rich product evaluations with species concentration and 
temperature measurements showed the attainment of a relatively uniform 
product composition. 

Phase 3 : Test Matrix and Data Grid Design 

The reacting experiment was designed as a continuation of the previous 
non-reacting experiments performed by Hatch, et al. (1996), Rroll, et al. (1996), 
and Sowa, et al. (1994). The experimental conditions were designed to simulate 
the conditions tested by Sowa, et al. where the momentum-flux ratio J was 40 
and the mass-flow ratio was 2.5. These momentum-flux and mass-flow ratios 
are representative of the design values for the proposed HSCT aircraft engines. 

Six modules with 8, 9, 10, 12, 14, and 18 orifices were used for the 
experiment. The round hole geometry was chosen as the baseline 
configuration. The choice of the 9-hole configuration was based on the non- 
reacting optimization results of Sowa, et al. which determined that the 9-hole 
case produced the best mixing at a momentum-flux ratio J of 40. The 8- and 10- 
hole cases were chosen to bracket the 9-hole case. In the course of testing, it 
was discovered that the optimum penetration was not obtained with the 9-hole 
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case, and successive tests were performed with the 12-, 14-, and 18-hole 
modules in order to encompass the overpenetrating to underpenetrating cases. 
The elapsed time to traverse to the specified grid point, to wait for the 

» 

readings to stabilize, and to obtain a datum point lasted between 1.5 to 2 

* 

minutes. The time constraint of operating the experiment in one continuous 
run per module necessitated the use of a coarse data grid across two orifices at 
six axial locations. Temperature and species concentrations were obtained to 
evaluate reaction and mixing occurring in the quick-mix module. 

Phase 4 : Analysis and Evaluation of Data 

The temperature and species concentration measurements obtained were 
processed and graphed to give a pictorial indication of the jet penetration and 
its effect on mixing and reaction. Recommendations for future tests and 
improvements were formulated based on the conclusions and experiences 
gained from the baseline tests and from bringing the reacting experiment to 


fruition. 


CHAPTER 4 


EXPERIMENT 

4.1 Facility 

The atmospheric model RQL combustion facility pictured in Figure 4.1 
was designed for fuel-rich reacting flow capability. Air and fuel flows to the 
up-fired facility were regulated through a flow panel. The experiment 
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Figure 4.1 Reacting Experiment Facility 
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consisted of a fuel-air premixing section, a refractory-lined stainless steel can 
which supported rich combustion, and a plenum-fed quick-mixing section that 
exhausted into a fume hood. Data were acquired via an intrusive probe that 

♦ 

was kept stationary while the entire test stand was traversed. 

* 

4.1.1 Air and Fuel Flow Supply 

Air and fuel flow regulation was accomplished through a centralized 
flow panel network shown schematically in Figure 4.2. A flow circuit 
previously used by Hatch, et al. (1996) and Kroll, et al. (1996) supplied the air 
flows to the facility. Dried and filtered air was supplied by an on-site air 
compressor factory. 

The main air flow to the rich combustion section and the air flow to the 
four independent jet air pathways were metered by sonic venturis. An 
additional rotameter was used to monitor fuel flow rates to the combustor. The 
fuel rotameter and sonic venturis were calibrated with a Laminar Flow 
Element. 

The choice of fuel for this experiment was based on operational and 
chemical considerations. The use of a gaseous fuel was desired in order to 
eliminate the complexities associated with liquid fuel atomization. Natural 
gas, though available to the experiment in ample supply, was not utilized 
because its highly refractory chemistry and rich flammability limit were not 
able to sustain a stable source of combustion products at an equivalence ratio <|> 
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above 1.3 in preliminary tests. The possible impurity and daily variation of 
natural gas composition also did not ensure a constant fuel composition to the 
experiment. 
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Figure 4.2 Schematic of Flow Panel (Adapted from Hatch, et al, 1996) 
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Propane (C 3 Hg) was chosen because its pyrolysis and thermal 
decomposition kinetics simulate the properties of jet fuel. Propane also 

possesses a rich flammability limit above an equivalence ratio <j> of 2.0, thus 
enabling the attainment of the rich equivalence ratio desired for the 

* 

experiment. The fuel was supplied by 94.6 L (25 gal.) tanks of liquefied 
propane. 

4.1.2 Rich Product Generation 

The challenge in designing the rich product generator for the experiment 
was to produce a consistent and uniform effluent of non-swirling, rich products 
into a cylindrical mixing module. Hardware durability and safety issues 
applied additional constraints on the design. 

The final design of the rich product generator incorporated in the 
experiment is shown in Figure 4.3. Propane and air at room temperature 
(25°C) were mixed along a 4.3 m length of a 5 cm diameter pipe before the 
ignition region. The ignition source was provided by an industrial spark plug 
placed in the center of a quarl section. The quarl provided a 3.8 cm contraction 
to prevent the backflow of combustion products. 

The pipe and combustor sections upstream and downstream of the quarl 
were cast with an aluminum oxide (A1 2 0 3 ) refractory material. The refractory 
material was cast in the upstream pipe section to form an inner diameter of 3.8 
cm, and cast in the downstream combustor section to form an inner diameter of 

I 
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8 cm to match that of the quick-mix module. The wall thickness of the cast 


ceramic compound was approximately 1.3 cm. The refractory material, which 


is rated up to 1870°C, insulates the reaction and prolongs the life of the stainless 

» 

steel combustor. 
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Figure 4.3 Rich Product Generator 
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The ignition procedure was formulated to eliminate the flashback 
concern associated with a premixed fuel-air system. The jet air was first 
supplied to cool the combustor and to prevent the reaction from exiting the 

♦ 

module through the orifices. The main air was then supplied to the system at a 
flow rate of 11.3 standard liters per minute (SLPM), which ensured a 10 m/s 
flow velocity through the quarl contraction. Assuming a turbulent flame speed 
of 2 m/ s, the velocity provided by the initial main air flow rate was sufficient 
to prevent flashback. The propane flow rate was increased and the spark plug 
switch was depressed until ignition occurred. The air flow rate was 
subsequently increased to a final setting of 22.7 SLPM, whereupon the fuel flow 
rate was increased to 1.51 SLPM. 

A recirculation region was needed to promote the stable continuous 
combustion of the fresh incoming fuel and air mixture. With the space 
constraints of the facility imposing a limit on the reaction residence time, the 
swirler offered the best solution to providing a compact recirculation region. 
Stable combustion was achieved through the use of a cast swirler with 45° 
vanes and holes dispersed around the outer circumference of the vanes. The 
design promoted mixing by impinging axial air flow through the holes with 

* 

the swirling air flow. 

Though the swirler aided in producing a stable reaction, a uniform plug 
flow was desired in order to avoid complications with data analysis and to 
provide a baseline case for future tests. An oxide-bonded silicon carbide 
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(OBSiC) ceramic foam matrix was used to dissipate the swirl imparted on the 
flow by the swirler. The ceramic is rated up to 1200°C and has been able to 
withstand thermal shock and stress loading due to periodic testing. The 7.62 
cm diameter, 2.54 cm-thick matrix was positioned such that it was five duct 
diameters downstream of the quarl section and one duct diameter upstream of 
the quick-mix module. The porosity of the foam, rated at 4 pores per cm (10 
pores per inch), was sufficient in allowing the required flow rate through the 
2.54 cm thickness at a negligible 0.3% pressure drop in the system. 

Profile measurements obtained across diameters of the rich zone showed 
that relatively uniform product concentrations and temperatures were 
achieved. The structural integrities of the ceramic foam and of the combustor 
refractory lining and parts were assessed after the apparatus was subjected to 
continuous testing lasting up to three hours per run. The apparatus has 
already undergone over 100 hours of total testing with no noticeable damage. 

In summary, the challenge of producing a rich product generator was 
met with the construction of a system that was safe to operate, exhibited part 
durability over continuous testing and cyclic loading, and produced a stable 
reaction and uniform rich effluent to the quick-mixing region. 

4.1.3 Jet Plenum Delivery 

The plenum feeding the jet air to the quick-mixing module was 
fabricated from aluminum pipe 15.2 cm in inner diameter to give a 3.4 cm 
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clearance between the plenum wall and quick-mix module. The air feeding the 
quick-mixing jets also served to cool the wall of the top two-thirds of the rich 
combustor length. Convective wall cooling is utilized in rich combustor wall 

% 

cooling (Zarzalis, et al, 1992) because conventional liner cooling schemes, 

* 

which introduce film air into the rich reaction, encourage the formation of hot 
stoichiometric pockets. 

Air for the quick-mixing jets was directed through flexible hoses to four 
air ports located near the bottom of the plenum. A high-temperature steel flow 
straightener 7.6 cm high with cell diameters of 0.95 cm was placed between the 
combustor and plenum walls to condition and promote an equal distribution of 
jet air entering the quick-mixing module. A pressure tap and a K-type 
chromel-alumel thermocouple monitored the pressure drop across the quick- 
mixing module and the temperature of the jet air, respectively. Optical access 
into the plenum was provided by two Pyrex windows situated 180° apart. 

4.1.4 Quick-Mixing Modules 

The modular quick-mixing section allowed for testing different jet 
orifice configurations by interchanging cylindrical quartz modules. The 

t 

modules were held in place by a sealing mechanism which compressed the 

* 

module against the combustor. Ceramic fiber paper consisting of a blend of 
alumina and silica that was rated up to 1260°C under continuous usage served 
as the gasket material between the quartz module and stainless steel surfaces. 
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The 80 mm inner diameter, 85 mm outer diameter modules were 280 
mm long, with the orifices equally spaced along the circumference of the 
tubing (Figure 4.4). The six geometries tested were the 8-, 9-, 10-, 12-, 14-, and 
18-round orifice configurations. 
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Figure 4.4 Quick-Mix Module Dimensions 


The total area for each orifice configuration was 1237 mm 2 and was 
based on a design momentum-flux ratio of 40 and discharge coefficient of 0.65. 

The calculation of the jet orifice area appears in Appendix A. 

* 

4.2 Data Acquisition 

Data were acquired with a stationary probe while the experimental 
apparatus was moved through the x-, y-, and z- traverses to the desired data 
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point. A digital encoder (Mitutoyo) monitored the position of the probe tip 
with respect to the center of the quick-mix module. A sector of data was 
obtained at each of six planes situated throughout the length of the module. 
Measurements of temperature and species concentration were obtained. 

4.2.1 Probe Design 

A double-jacketed water-cooled stainless steel probe 30 cm in length was 
used to extract gas samples from the quick-mixing section (Figure 4.5). The 



Figure 4.5 Probe Design 


38 


probe measured 8 mm in outer diameter and tapered down to a 3.2 mm tip. A 
45-degree bend was made one inch from the tip. The probe design was 
influenced by the research of Sowa, et al. (1994), who found that a 

* 

thermocouple probe with a 45-degree angled tip was the best design for 
acquiring temperature data that biased the mainstream and incoming jet flow 
equally in the orifice region. The plane of the angled probe tip was positioned 
such that the tip pointed toward the center of the sector wall. 

4.2.2 Data Grid and Planes 

Six cross-sectional planes of data were obtained per module as depicted 
in Figure 4.6. With z referring to the axial distance, R defined as the module 
radius, and d defined as the orifice diameter, the planes were situated, with the 
origin z = 0 set at the leading edge of the orifice, at positions 

(1) one module radius upstream (z/R = -1), 

(2) at the orifice leading edge (z/R = 0), 

(3) one-half the orifice axial height (z/R = (d/2)/R), 

(4) at the orifice axial height (z/R = d/R), 

(5) one module radius downstream (z/R = 1), 

* 

(6) two module radii downstream (z/R = 2). 

Assuming flow symmetry, a sector of data was obtained to represent the entire 
plane (Figure 4.7). Each plane of data consisted of 16 points distributed across 
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a sector that encompassed two orifices. One point was located at the center 
while the rest of the points were positioned along the arc lengths of three radii 
at 12.7 mm, 25.4 mm, and 38.1 mm. Five points were positioned along each arc 
length such that their positions formed radial-axial planes that bisected either 
the orifices or the wall region between the orifices. Of the resulting radial-axial 
planes, two were aligned with the orifices and three were aligned with the 
wall. 

A 90° sector was used for the 8-hole module, an 80° sector for the 9-hole 
module, a 72° sector for the 10-hole module, a 60° sector for the 12-hole 


module, a 51.4° sector for the 14-hole module, and a 40° sector for the 18-hole 
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module. The data grid density increased with the number of orifices in the 
module because the size of the sectors decreased. 




Figure 4.7 Data Grid Sector over Two Orifices 


4.2.3 Temperature Measurements 

A B-type platinum-rhodium thermocouple was constructed from a set of 
bare 30% versus 6% Rhodium wires 0.254 mm in diameter. The range of the B- 
type thermocouple, which falls between 0° and 1700°C, provided a suitable 
range for the reaction temperatures that were measured. The error associated 
with the thermocouple wire was 0.5°C for temperatures above 800°C. 

Both wires, with one ensconced in a cloth sheath, were inserted into 
plastic tubing that was shrunk though heat-treatment. The thermocouple end 
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that would be inserted into the reacting flow had its wires threaded separately 
through two channels of a 2.5 cm length of ceramic tubing. The ends of the 
wire exiting the ceramic were joined with a spot weld. 

The thermocouple had a thickness that was less than 3 mm to facilitate 
its threading through the probe. The unjoined ends of the thermocouple 
exiting from the probe were connected to an analog-to-digital screw terminal 
panel attached to a personal computer. An ice water bath served as the 
reference point for the thermocouple. Figure 4.8 depicts a schematic of the 
temperature acquisition set-up. 
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Figure 4.8 Temperature Data Acquisition Set-Up 


A data acquisition program (Omega Engineering, Inc.) was modified to 
read the voltage signal from a B-type thermocouple. The program sampling 
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rate was set at approximately 100 readings per 15 seconds. Voltage 
measurements were obtained at a point as soon as the readings stabilized after 
a period of approximately 45 seconds. The fluctuation in temperature 
measurement was highest in the jet orifice region, where fluctuations reached 
at most 100°C, or 10% of the mean measured value. 

The readings acquired by the program were compared to readings 
obtained from a digital meter (Analogic AN2402) at several points in the rich 
section of the combustor. The difference between the readings was determined 
to be less than 5%. 

4.2.4 Species Concentration Measurements 

Species concentrations by volume of CO, CO z , 0 2 , unburned HC, and 
NO x were obtained by drawing gas samples from the points in the flowfield 
planes. The samples were drawn by a vacuum pump from the probe through a 
15 m heated line connected to the emission analyzers (Figure 4.9). The heated 
line was maintained above 100°C to prevent water drop out in the line and 
subsequently the possibility of water-soluble NO x dissolving in the condensate. 

The analyzers (Horiba Instruments, Inc.) utilized non-dispersed infrared 
(NDIR), paramagnetic, flame ionization detection (FID), and 
chemiluminescence techniques to measure both C0 2 and CO, 0 2 , unburned HC, 
and NO x , respectively. Appendix B describes further the principles behind the 
detection methods employed by the emission analyzers. 
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Figure 4.9 Emission Analyzer Routing (Courtesy of M. Miyasato) 


Prior to testing, the analyzers were set to zero by flowing nitrogen gas 
(or air for the FID analyzer only) through the system. The analyzers were then 
calibrated with their respective span gases. 

A measurement was obtained at a point after the readings on the console 
stabilized (approximately 45 seconds). Data were acquired with an acquisition 
program that read 100 samples in 20 seconds and returned an averaged 
quantity. The uncertainty in the analyzer species concentration measurement 
was 1% of the full scale reading. 

4.3 Experimental Conditions 

The operating conditions under which the tests were run are noted in 
Table 4.1. The actual momentum-flux ratio was higher than expected because 
the orifice area was based on a rich zone adiabatic flame temperature of 1800K 
(see Appendix A) while the measured temperature was 20% lower. The jet 
temperature had also been underestimated, but the effect of the jet temperature 
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on the orifice area calculation was not as great. In addition, the actual orifice 
discharge coefficient C d of 0.73 was greater than the estimate of 0.65 used to 
design the orifice area. Appendix C shows the calculations made in deriving 
the values listed in Table 4.1. 
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Table 4.1 Operating Conditions 


Parameter 

Value 

P(kPa) 

101 

rich equivalence ratio <j) 

1.66 

overall (j) 

0.45 

■WK) 

1500 

T ( „ (K) 

480 

(m/s) 

18 

momentum-flux ratio J 

57 

mass ratio MR 

2.5 

density ratio DR 

3.3 

velocity ratio VR 

4.2 

discharge coefficient C d 

0.73 
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CHAPTER 5 


RESULTS AND DISCUSSION 


5.1 Overview 

Temperature and species concentration measurements were obtained for 
8-, 9-, 10-, 12-, 14- and 18-hole configurations of equal orifice area at a set 
momentum-flux ratio of 57 and mass-flow ratio of 2.5. The results are 
presented in three different formats: (1) histograms to depict the distribution of 
the raw data, (2) contour plots of radial-axial sections to show the axial 
evolution of the flow, and (3) contour plots of sectors to show planar symmetry 
and the extent of mixing and reacting processes. 

Table 5.1 summarizes the z/R values for each of the six planes measured 
per module. Planes 3 and 4 are situated at the orifice mid-height and height 
levels. Differing z/R values for each configuration occur because the orifice 
area is kept constant to keep the jet-to-mainstream momentum-flux ratio set. 


Table 5.1 Normalized Axial Length z/R per Plane per Module 


Plane 

8 Holes 

9 Holes 

z/R 
10 Holes 

Values 
12 Holes 

14 Holes 

18 Holes 

1 

-1 

-1 

-1 

-1 

-1 

-1 

2 

0 

0 

0 

0 

0 

0 

3 

0.175 

0.165 

0.157 

0.143 

0.133 

0.117 

4 

0.350 

0.331 

0.314 

0.286 

0.265 

0.234 

5 

1 

1 

1 

1 

1 

1 

6 

2 

2 

2 

2 

2 

2 
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5.2 Reacting Flow Field Description: 8-Hole Module 

The description of the reacting flow field of the 8-hole case is discussed 
first to gain a sense of jet mixing and reaction in a rich reacting crossflow. 

5.2.1 Temperature Profiles 

Figure 5.1 depicts histograms of temperature measurements obtained for 
the 8-hole configuration. The histograms provide information on the 
temperature distribution per plane at 16 grid points (recall Figure 4.7). Point 1 
lies in the center of the module. Points 2-6 lie along the arc length at one-third 
of the module radius (hereafter denoted as RJ, Points 7-11 lie at two-thirds of 
the module radius (denoted as RJ, and Points 12-16 lie near the module radius 
R. Points 3, 5, 8, 10, 13, and 15 are aligned with the orifices. 

Plane 1 (z/R= -1), which is representative of the rich product entering 
the quick-mixing module, shows a relatively uniform temperature distribution. 
Temperatures near 1500K were measured at points at the center and on the R, 
and R, arcs. Slightly lower temperatures ranging from 1300 to 1450K were 
obtained at the wall. The lower temperatures are attributed to the convective 
cooling of the outer module wall by the plenum air. 

At the orifice leading edge (Plane 2, z/R=0), the temperatures still hover 
near 1500K with the exception of deviations at Points 13 and 15. The lower 
temperatures occur since the points are in the near field of the jet entrance. 
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Figure 5.1 Temperature Distribution for the 8-Hole Module 
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By the orifice half-plane (Plane 3, z/R=0.175), the temperatures at a 
majority of points remain at 1500K. However, jet entrance is clearly seen with 
lower temperatures near 400K occurring at Points 13 and 15. Jet penetration to 
the second radius R, is marked by lower temperatures of 810K and 1300K 
occurring respectively at Points 8 and 10. 

The unequal temperatures measured at Points 8 and 10 should not occur 
if jet and crossflow symmetries are assumed. The symmetry assumptions 
appear valid because of the flow conditioning methods utilized in the 
experiment: a flow straightener installed in a long plenum chamber was used 
to ensure an equal distribution of jet flow into the module while a porous 
ceramic matrix was used to condition and promote a uniform reacting 
crossflow. 

To investigate the cause of the asymmetry, pitot tube measurements 
were obtained under non-reacting conditions at the entrance of both orifices 
and along two diameters across the inlet duct to the quick-mixing section (at the 
z/R= -1 plane). The measurements at the orifices showed a 3% to 5% velocity 
variation which suggests that the jet flow entering the two orifices is nearly 
equal. Measurements of the crossflow showed a near-uniform velocity profile 
(within a 15% range). Because Points 8 and 10 at Plane 3 lie in the vicinity of the 
jet-crossflow interface, the unequal temperature measurements could be 
attributed to a fluctuation created by fluid interaction and reaction processes. 
Time-resolved measurements are needed to verify this hypothesis. 


49 


At the orifice trailing edge (Plane 4, z/R=0.350), all of the jet fluid has 
entered the module. The jets penetrate to the second-third mixer radius R, but 
have not yet fully dispersed and mixed with the neighboring node points. Note 
that the central core temperatures are still maintained near 1500K. 

Substantial reaction and mixing occurring between the plane at the 
orifice trailing edge (Plane 4) and the plane one duct radius downstream of the 
orifice leading edge (Plane 5, z/R=l) is suggested by major differences between 
the two histogram profiles. At Plane 5, mixing and reacting processes have 
produced the following evenly distributed, stratified bands of temperature: a 
1000K band at the center and first mixer radius R,, a 1400K band at the second 
mixer radius R^ and a 1750K band near the module radius R. The 1750K 
temperatures in the outer band R are higher than the initial rich effluent 
temperature of 1500K, which suggests the occurrence of chemical reaction in 
that region. 

Jet fluid penetration to the central core is inferred by the temperatures 
observed in the central core which are overall lower than the 1500K 
temperature of the entering rich effluent. The 1000K temperature band at the 
center and at R, are higher than the initial 400K temperature of the jet fluid, 

i 

which suggests the occurrence of either reaction, mixing, or both processes. 
However, the extent to which the higher temperatures are attributed to reaction 
rather than to mixing cannot be determined without measuring conserved 
scalar quantities to derive the mixture fraction field. 
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The histogram profile obtained two duct radii downstream of the orifice 
leading edge (Plane 6, z/R=2) is nearly identical to that obtained at Plane 5. 
The similar profiles, which show a lack of extensive mixing and reacting 
activity between the two planes, indicate that the bulk of the reaction and 
mixing profiles appears to have been completed within one duct radius of the 
jet entrance. 

The contour plots in Figures 5.2a and 5.2b offer a spatial view of the 
temperature distributions. The radial-axial and sector cross-sectional plots 
respectively depict the longitudinal and planar evolution of the flow field. 

The radial-axial plots (e.g. in Figure 5.2a) depict two cross-sections: one 
which is aligned with the wall midpoint between the orifices (a "midplane" 
cross-section) and one which is aligned with the orifice (a "centerplane" cross- 
section). The midplane cross-section is an average of the three lines of data 
aligned with the wall while the centerplane cross-section is an average of the 
two lines of data aligned with the hole. The data are plotted on axial versus 
radial length scales z and r normalized with respect to the mixer radius R. 
Sector plots such as those seen in Figure 5.2b are an unaveraged depiction of 
the interpolated 16 points of data obtained per sector. 

A nearly uniform temperature distribution in the rich flow, with the 
exception of the small band of lower temperatures occurring at the wall, enters 
the quick-mixing regime. Jet introduction into the flow is seen at the orifice 
mid-plane (Plane 3, z/R=0.175), with the full introduction of jets occurring by 
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the orifice trailing edge (Plane 4, z/R=0.350). The longitudinal centerplane plot 
(Figure 5.2a) depicts jet penetration toward the centerline within one duct 
radius (Plane 5, z/R=l). The midplane cross-section shows relative symmetry 
with the centerplane cross-section in the flow downstream of the z/R=l plane. 

Directly downstream of the orifices, temperatures are lower than the 
initial mainflow temperature as the cooler jet fluid entrains and mixes with the 
hot crossflow fluid. At one duct radius downstream (z/R=l) the jets penetrate 
toward the center and displace the rich crossflow to the wall. The displaced 
fluid undergoes an increase in temperature as some of the jet fluid mixes with 
the crossflow to initiate the CO oxidation reaction. As a result, a band of 
temperatures higher than the initial rich effluent is formed. 

The sector cross-sections in Figure 5.2b show the occurrence of 
temperature stratification beginning at Plane 5. At two duct radial lengths 
(Plane 6, z/R=2) the stratification becomes more distinct but does not change 
drastically in distribution. 

5.2.2 Species Concentration Profiles 

Figures 5.3a and b show the histogram distributions of 0 2 , C0 2 , CO, HC, 

V 

and NO x concentrations measured for the 8-hole module. The concentrations 
are measured in terms of either % or parts per million (ppm) by volume. The 
distributions at Plane 1 reveal uniform rich zone concentrations of 0% 0 2 , 5.2% 

C0 2 , 12% CO, and 1.5 ppm NO x . The HC concentrations, while overall 
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relatively uniform, show more variability. 1% to 1.6% concentrations occur in 
the central core while higher concentrations up to 2% are present in the wall 
region. HC species are comprised of products of incomplete combustion (PICs) 
that have not oxidized because of either a lack of oxidant or temperature to 
propagate the reaction. The latter explanation applies in this case, where the 
lower wall temperatures seen in Figures 5.1, 5.2a, and 5.2b promote the higher 
HC levels near the wall region. 

Mass balances of carbon, hydrogen, oxygen, and nitrogen were 
performed on the rich products measured at Plane 1 to determine the integrity 
of the measurements. The reaction of propane (C 3 H 8 ) and air was assumed to 
yield products of CO, C0 2 , unburned C 3 H 8 , 0 2 , NO, N 2 , and H,0. The fuel and 
air flow rates are known; CO, C0 2 , unburned C 3 H 8 , 0 2 , and NO are measured; 
and N 2 is assumed to make up the rest of the exhaust gas volume. The only 
unknown quantity in the reaction is the mole fraction of H,0, which has been 
dropped out from the sample by a condenser before the sample enters the 
emission analyzers. As a result, the mole fraction of water is calculated from a 
separate mass balance performed on each element and is compared for 
equivalency. 

* 

The resulting analysis (see Appendix D) shows the water mole fraction 
for all elements except carbon to be within the same range. However, a 
comparison between the measured species concentrations obtained in the rich 
product region with output from the NASA chemical equilibrium code shows 
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similar results: under room temperature and atmospheric pressure the NASA 
equilibrium code predicts that the combustion of propane at an equivalence 
ratio of 1.7 yields 0% 0 2 , 4.1% C0 2 , 12.3% CO, and 1.3 ppm NO x . Another major 
constituent that the equilibrium code predicts but that was not measured in the 
experiment is H,, for which a 9.4% concentration was predicted. The 
neglection of PL, in the mass balance equation may have greatly affected the 
closure of the equation. 

Inferences related to jet penetration and mixing can be formed from the 
axial distribution of species concentrations. The 0 2 charts in Figure 5.3a, for 
example, show the evolution of jet penetration from the leading edge of the 
orifices (Plane 2). Farther downstream at the orifice mid-height plane (Plane 3), 
the jets fill the outer mixer radius R to near-21% concentrations at Points 13 and 
15, which are aligned with the orifices. The jet fluid then migrates toward the 
second radial band Rj (Plane 4), as seen by 21% O, concentrations at Points 8 
and 10, and begins to disperse, as noted by the appearance of oxygen at 
neighboring points 7, 9, and 11. At one and two mixer radii downstream of the 
orifice leading edge (Planes 5 and 6) the jet fluid penetrates to the central core. 

At the orifice leading edge and mid-orifice planes (Planes 2 and 3) the 
points at which 0 2 levels are measured (Points 13 and 15 in Plane 2, Points 8, 10, 
13, and 15 in Plane 3) correspond to the decrease or increase in the other species 
concentrations measured at the same points. The decrease in C0 2 levels at 
Points 13 and 15, and at Points 8 and 10 at Plane 3 is attributed to jet dilution 
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(see Figure 5.3a). The slight increase of C0 2 concentrations at other points in 
the plane reflect reaction. Similarly, the dilution of CO (Figure 5.3a) and HC 
(Figure 5.3b) at Points 13 and 15 correlates to the presence of 0 2 concentrations 
at the same points in Plane 2. 

The occurrence of reactions that consume CO and HC is corroborated by 
a comparison between the overall concentrations measured at Planes 2, 3, and 4. 
Plane 3 shows a decimation of HC concentrations that had been present in 
abundance in Plane 2. The CO concentrations are also decreased, though not as 
rapidly as HC, because of the slower CO reaction. Between Planes 3 and 4, 
additional CO has reacted while no further reaction of HC occurs. 

The presence of reaction in Planes 2 and . 3 is verified by the NO x 
histograms (Figure 5.3b). At Plane 2, Points 13 and 15 show slightly higher NO x 
concentrations, a sign that reaction processes occurring at the jet-crossflow 
interface are increasing temperatures to levels that promote the formation of 
NO x from the nitrogen in the air. The reaction continues in Plane 3 where 
points at the second radius R, register higher levels of NO x . Between Planes 3 
and 4, however, jet displacement processes dominate as the body of jet fluid is 
introduced into the crossflow. No further evidence of substantial NO x 

¥ 

production is seen; the NO x produced in the second radial region migrates to 
the first radius R and the centerline. 

At planes within the vicinity of the orifices (Planes 3 and 4), all of the 
measured species show asymmetric concentrations in the second radial band R 
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where the points aligned with the orifices (Points 8 and 10) and the points 
aligned with the wall (Points 7 , 9, and 11) should show equal concentrations 
within each group. The asymmetric species concentrations measured in the 
second radial band may be attributed to fluctuations caused by reacting and 
mixing processes occurring at the jet-crossflow interface, as explained for the 
unequal temperature distribution similarly noted in section 5.2.1. 

The bulk of the reaction takes place in the zone downstream of the holes 
(between Planes 4-5) where the jets penetrate toward the center and begin to 
disperse throughout the crossflow. In this region of reaction between Planes 4 
and Plane 5, C0 2 and NO x concentrations rise, particularly in the wall region, 
while CO and HC concentrations decrease significantly. Most of the reaction 
terminates by the z/R=l plane (Plane 5) since the magnitude of the profiles 
measured at that plane and at z/R=2 (Plane 6) do not change substantially. 

The contour plots in Figures 5.4a and 5.4b give an alternative 
presentation of the species concentration profiles. The radial-axial profiles 
(Figure 5.4a) illustrate the evolution of jet penetration and flow field 
characteristics in the module while the sector plots (Figure 5.4b) offer another 
view which indicate species distribution symmetry. 

From Figures 5.4a and 5.4b, the species entering the module are 
generally uniform up to the plane of jet entry. Evidence of early jet penetration 
is apparent at the orifice leading edge (z/R=0). Full jet penetration is apparent 
for all species by the orifice trailing edge (z/R=0.350). 
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Jet penetration as seen in the 0 2 profiles is represented by a high 
concentration of jet fluid in the orifice region that enters the crossflow nearly 
intact before dispersing throughout the mixer radius (see Figure 5.4a). The jet 
appears to penetrate toward the center within one mixer radius length (i.e., at 
Plane 5, z/R=l). 

The centerplane profiles show a transition region in the jet wake formed 
by the initial reaction between the jets and the entrained rich crossflow. 
Downstream of the z/R=l plane the penetrating jets displace the rich reacting 
fluid toward the wall, as evidenced by the pocket of high CO concentration at 
the wall. The small pocket of CO subsequently disappears as CO reacts with 
the jet to form C0 2 . This source of C0 2 increases concentrations at the wall 
from the z/R=l to z/R=2 plane. The higher NO x values also occurring at the 
wall downstream of the z/R=l plane correspond to the higher temperatures 
(recall Figure 5.2a) that are associated with the CO oxidation reaction. 

In addition to showing the extent of reaction and of mixing occurring 
between the jets and the rich mainstream, the sector plots in Figure 5.4b give an 
indication of jet and flow field symmetry. Although at the orifice midplane 
(z/R=0.175) the first jet appears to dominate over the second jet, by the orifice 
trailing edge (z/R=0.350) the jet flows are essentially symmetrical. At the 
z/R=l plane, the flow field of stratified concentration bands is radically 
different from the previous plane which shows discrete jets. At the z/R=2 
plane the stratified bands become more well formed but the overall flow field 
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does not change drastically. This observation as well as the virtual 
disappearance of CO and HC between the z/R=l and z/R=2 planes show that 
the jets have almost entirely reacted with the rich crossflow within one duct 
radius of jet entry for the 8-hole configuration. 

5.3 Comparison of 8-, 9-, 10-, 12-, 14 -, and 18-Hole Modules 

The mixing and reaction performance of the six orifice configurations is 
assessed with a comparison of radial-axial and sector contour plots for each 
module. The objective of comparing the six mixing modules is to determine a 
configuration that leads to optimal jet trajectory penetration and hence optimal 
mixing and reaction within one duct radius of the jet entrance (i.e., at z/R=l). 
In combustor design, the attainment of a short combustor length is desirable to 
maintain the compactness of the engine. For the RQL configuration it is thus 
preferable to attain complete mixing in a minimal length. The one duct radius 
limitation has been used as an arbitrary reference plane of comparison by 
Hatch, et al. (1996), Kroll, et al. (1996), and Sowa, et al. (1994). 

For the module performance comparison, only the contour plots are 
presented here. Stacked histogram plots of the temperature and species 
concentration profiles for the other module configurations (9-, 10-, 12-, 14-, 18- 
hole) may be found in Appendix E. 
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5.3.1 Temperature Profiles 

Figure 5.5 shows the axial evolution of temperature profiles at the center- 
and midplane sections of all the cases tested. All six cases show rich inlet 
product temperatures of 1500K as well as cooler wall temperatures between 
z/R= -1 and z/R=0. The modules share other similar flow features such as a 
400K temperature occurring where the jet enters the orifices and a band of hot 
fluid at the wall downstream of z/R=l. 

The jet trajectory, defined by the approximate locus of points showing 
minimal temperature with respect to distance, is represented on the centerplane 
plots by a red line extending from the orifice mid-height to the z/R=l plane. 
Beyond the z/R=l plane, the disintegrating jet structure makes the 
determination of the trajectory difficult. 

As the number of holes is increased, the jet penetration as indicated by 
the trajectory decreases and the jet flow disperses into the crossflow rather than 
accumulates in the central core. If the path of the trajectory is extended past the 
z/R=l plane, the 8- through 12-hole modules produce jet impingement at the 
centerline and varying degrees of jet overpenetration. Jet trajectory inclination 
toward the center causes more fluid to concentrate in the center of the mixer at 
the z/R=2 plane, which is seen in the 8- to 12-hole cases by the cooler fluid 
temperatures there. The 14-hole case shows the attainment of a near-uniform 
temperature distribution between 1300 and 1600K (orange bands) by the z/R=2 
plane whereas the 18-hole produces a hot core throughout the length of the 
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module which suggests jet underpenetration. Hence, the uniform temperature 
distribution produced by the 14-hole case shows that the configuration is 
optimal for this experiment. In this case, the jet trajectory reaches the half- 
radius of the mixer by the z/R=l plane. 

The sector profiles in Figure 5.6 show a stratification of fluid 
temperatures occurring by the z/R=l plane which becomes more concentrically 
formed by the z/R=2 plane. In each case, there is an increase in temperature 
occurring between the z/R=l and 2 planes, which suggests the occurrence of 
further reaction beyond one duct radius of jet entry. The increase is not as 
apparent in the 8-, 9-, and 10-hole cases as it is in the 12-, 14-, and 18-hole cases, 
which indicates more reaction occurring in the latter three configurations. 

In the 8-hole case, the presence of fluid between 800K and 900K occurs in 
the center of the z/R=2 plane whereas the center temperatures at this plane for 
the 9-, 10-, and 12-hole modules are between 900K and HOOK. The 18-hole case 
shows a reverse trend where the center of the z/R=2 plane possesses hotter 
temperatures ranging between 1600 and 1750K and the outer annular region 
contains cooler fluid temperatures between 1150 and 1450K. The 14-hole case 
exhibits the most evenly dispersed bands of temperatures at the z/R=2 plane. 
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5.3.2 Species Concentration Profiles 

The effect of the number of holes on jet trajectory penetration is better 
illustrated in the radial-axial 0 2 distributions of Figure 5.7. In the centerplane 
plots of all cases, discontinuities in color bands (yellow to blue bands) 
emanating from the jet orifice to the z/R=2 plane occur. The re-emergence of 
bands of high species concentration suggests the implausible production of 
additional 0 2 in the reacting system. The lack of data measured between the 
orifice trailing edge and the z/R=l planes, and between the z/R=l and z/R=2 
planes probably accounts for the discontinuous bands of 0 2 concentration 
produced by interpolation. 

The jet trajectory, defined as the locus of points showing maximal 0 2 
concentration as a function of distance, is drawn on the centerplane plots from 
the orifice midplane to the z/R=l plane. The trajectories based on 0 2 
concentration nearly coincide with those based on temperature (see Figure 5.5). 

High penetration of the jet trajectory causing fluid impingement at the 
centerline occurs in the 8-hole case and results in a large concentration of 0 2 in 
the central core. This condition suggests that the 8-hole module is an 
overpenetrating case which is undesirable because the oxygen tends to 
accumulate in the center rather than disperse, mix, and react with the crossflow. 
Overpenetration also leads to less reaction since the accumulation of jet fluid in 
the center decreases the area of the jet-crossflow interface, or the promotion of 
reactant interaction. The 9-, 10-, and 12-hole cases also exhibit some degree of 
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overpenetration as the jet trajectory appears to intersect with the centerline by 
the z/R=2 plane. 

The 14- and 18-hole cases produce a jet trajectory that does not penetrate 
to the center by the z/R=2 plane. Without jet impingement, a more lateral 
spreading of jet fluid, as evidenced by the narrower distribution of oxygen 
concentration values across the mixer radius, is achieved. The midplane 0 2 
profiles show that the 14- and 18-hole modules exhibit a smaller band of fluid 
devoid of 0 2 at the wall (indicated in dark red and orange) between the z/R=l 
and z/R=2 planes. The 18-hole case in fact shows no red bands between 0-4% 
in the jet wake, unlike the 14-hole case. However, the jet underpenetration 
property of the 18-hole case manifests itself with the occurrence of lower 0 2 
concentration bands in the core region. Like jet overpenetration, jet 
underpenetration decreases the maximal jet-mainflow surface area of reaction 
because a portion of the jet fluid is bounded by the wall of the module. Jet 
underpenetration allows the rich reaction products to exit the module without 
completing the combustion process, which is undesirable. In the 14-hole case, a 
larger jet surface area exposed to the crossflow accelerates jet dispersion and 
reaction such that more of the rich mainstream flow reacts with rather than 
bypasses the jets. 

For cylindrical crossflow geometries, several investigations have 
determined a jet penetration depth that leads to better mixing. In a numerical 
study performed by Talpallikar, et al. (1992), results suggest that optimal mixing 
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occurs when the jet penetrates to the mid-radius. Kroll, et al. (1993) infers from 
experimental results that optimal mixing occurs when the jet penetrates to the 
radius that divides the mixer into an equal core and annular area, or at a radial 
distance 30% from the wall. For the cases tested in this experiment, the jet 
trajectory as determined by both oxygen and temperature measurements 
intersects the z/R=l plane beyond the radial midpoint from the wall for the 8-, 
9-, and 10-hole cases; appears to approach the mid-radius for the 12- and 14- 
hole cases; and lies at a point before the mid-radius and toward the wall for the 
18-hole case. The Talpallikar, et al. criteria support either the 12- or 14-hole case 
as producing the jet penetration that promotes the best mixing out of the 
configurations tested while Kroll, et al. (1993) data supports the 14- or 18-hole 
case as being the optimal configuration. 

The number of holes leading to optimal jet penetration can be predicted 
by an empirical relation developed for non-reacting jets injected into a subsonic 
cylindrical crossflow. The relation from Holdeman (1993) states that the 
appropriate number of holes n that will lead to optimal penetration may be 
determined by the following equation (5.1): 


n = 


tcV2J 

C 


(5.1) 


where J is the jet-to-mainstream momentum-flux ratio, and C is a constant 
whose value is 2.5 if an optimum, single row injection is desired. Note that this 
relation includes the assumption that the "optimum" spacing for a rectangular 
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duct applies at the radius that divides the can into equal area can and annular 
sections (Holdeman, 1993). In the reacting experiment where J=57, the equation 
yields an optimal configuration of 13.4 holes. This calculation corroborates the 
designation of the 14-hole case as the optimal reacting case tested, but also 
suggests that distributions of an even higher if not comparable uniformity may 
be attained for a 13-hole configuration. 

The corresponding cross-sectional sector plots for 0 2 (Figure 5.8) show 
the jets entering the crossflow symmetrically for each module. It is observed 
again from the 0 2 sector plots that increasing the number of holes lowers jet 
penetration. By the orifice trailing edge all of the jet fluid should have entered 
and been accounted for in Plane 4. The sector cross-sections show a larger bulk 
of jet fluid occurring per orifice in the 8-hole case and a successive decrease in 
bulk jet fluid per orifice as the number of orifices increases. This decrease in 
mass flow per orifice is attributed to the decrease in area per orifice, because the 
jet velocity through each orifice is constant for all six configurations. 
Decreasing the individual jet mass flow rate decreases the jet momentum, 
which consequently diminishes jet penetration into the crossflow. 

The stratified concentration distributions at z/R=l (Plane 5) differ 
substantially from Plane 4, which is situated at the orifice trailing edge. For the 
overpenetrating cases (8, 9, 10, 12 holes), reaction as seen by a decreased 0 2 
concentration occurs beyond the z/R=l plane but is not as substantial as the 
further reaction that occurs beyond this plane for the optimal 14-hole and the 
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underpenetrating 18-hole cases. Increasing the number of holes increases the 
total circumference, or the total available surface area of reaction between jet 
and rich reacting product, and this coupled with a near-optimal degree of 
penetration results in further reaction beyond z/R=l. The 18-hole case, though 
producing underpenetrating jets, does not produce a severe case that cause the 
jets to flow along the module wall and decrease the availability of the jet surface 
area that allows for further reaction. 

The experiment was designed to transition from a <|>=1.7, fuel-rich section 
(seen as Plane 1) to a (j)=0.45, fuel-lean section (seen downstream of the jet 
orifices as Planes 5 and 6). For propane combustion at an equivalence ratio 
<|)=0.45 at room temperature and pressure, the NASA equilibrium code predicts 
an 0 2 concentration of 11.1%. From the scale for the 0 2 concentration profiles, 
11.1% falls within the yellow-green color band, which spans from 10-12% 
concentration. The 0 2 sector plots in Figure 5.8 show that by the z/R=2 plane, 
the 14- and 18-hole configurations yield larger areas of the 10-12% color band. 

A comparison of the C0 2 concentration profiles for all six hole 
configurations (Figures 5.9 and 5.10) yields similar observations gained from 
the 0 2 species profiles in regard to the distribution of species concentrations. 
From the axial profiles (Figure 5.9) either the 14- or the 18-hole case appears to 
produce a more evenly dispersed C0 2 field downstream of the orifices. The 
sector profiles (Figure 5.10) offer a better viewpoint which shows the 14-hole 
case containing the least amount of gradient bands, or the most uniform profile 
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at the z/R=2 plane. The NASA equilibrium code prediction of a 5.47% C0 2 
concentration for <|)=0.45 falls within the yellow color band on the C0 2 plots. In 
the sector plots the 14-hole module exhibits a larger area with the yellow color 
band. 

Both the 14- and 18-hole cases appear to produce the highest overall 
amount of C0 2 distributed in the sector plane at the z/R=l and z/R=2 planes, 
which suggests that more complete reaction processes have occurred. In 
addition, both of these cases also show a higher degree of reaction occurring 
past the z/R=l plane that was noted earlier in the 0 2 sector plots. 

A sense of the extent of reaction can be gleaned from the radial-axial CO 
profiles (Figure 5.11). Between z/R=l and z/R=2, the presence of CO in the 
wake of the jet coupled with an increase of C0 2 in the same region suggests that 
the reaction of CO is a major contributor to the increased C0 2 at the wall (see 
Figure 5.9). As the number of orifices is increased, the pocket of rich CO-laden 
gases in the jet wake decreases in size and concentration until it is all nearly 
consumed by the 18-hole case. However, only the 18-hole case shows a higher 
CO concentration between 2.6 and 5.2% in the core by the final z/R=2 plane. 
The evidence of a rich product core in the 18-hole case supports its 
underpenetration designation. 
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Figure 5.12 shows the planar cross-sectional view of the CO distribution. 
For all of the modules a majority of CO has been consumed by the z/R=l plane. 
However, CO-rich pockets in the jet wake are observed in the same plane for 
the overpenetrating 8- through 10-hole cases, a CO-rich core is seen in both the 
optimal 14-hole and the underpenetrating 18-hole cases, and a balance of both 
is seen in the overpenetrating 12-hole case. Curiously, by the z/R=2 plane, only 
the 12-hole case exhibits a uniform, low CO-concentration band less than 1.3%. 
The NASA equilibrium code prediction of CO concentration from atmospheric 
propane combustion at <{>=0.45 is in the ppb range and can be considered 
negligible, which leaves only the 12-hole case with the achievement of a lower 
overall CO concentration by the z/R=2 plane. 

The HC contour plots (Figures 5.13 and 5.14) proffer little significant 
information. A comparison of HC profiles for each module in both axial and 
sector views shows a near-total consumption of HC which leads to a near-total 
concentration distribution of HC below 1% by the orifice trailing edge (Plane 4). 

A NO x profile comparison (Figures 5.15 and 5.16) yields the same general 
observations made in the other species profile comparisons. Similar to the C0 2 
species distribution, the bulk of NO x formation appears downstream of the 
orifices along the wall region for all cases, but also along the center region for 
the 14- and 18-hole cases. Among the six modules, the 14-hole case produces 
the most uniform NO x concentrations across the z/R=l and z/R=2 planes (see 
Figure 5.16). 
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The measured concentrations of NO x were found to be less than 40 ppm 
per point. It is emphasized, however, that the NO x concentrations measured in 
this experiment may not be indicative of thermal NO x behavior because the 
experiment is not run with air preheat, and peak temperatures are therefore 
suppressed. Any correlation made at this juncture between jet mixing and NO x 
production would not be reflective of situations occurring at actual engine 
operating conditions. 
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CHAPTER 6 


SUMMARY AND CONCLUSIONS 


6.1 Conclusions 

An experiment has been designed and successfully demonstrated to 
provide a test bed for the study of jet mixing in a rich reacting environment. In 
this initial demonstration, it was possible to determine the jet trajectory, as well 
as mixing and reacting processes, for six round hole configurations. 

For a set momentum-flux ratio J = 57 and mass-flow ratio MR = 2.5 it 
was found (under atmospheric conditions) that: 

• The data grid density and planar measurement distribution provide 
sufficient information to form general inferences and comparisons of 
mixing and reacting properties between various multiple-orifice 
configurations. 

• Profiles of temperature and species concentration can be used to an 
extent to indicate general zones where mixing and reacting processes 
occur. However, the separation of reaction from mixing cannot, at 
this juncture, be determined. A conserved scalar measurement of, for 
example, an inert gas tracer will enable a pure characterization of 
mixing in the reacting flow field. 
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• The jet wake is the site of further reaction downstream of the one 
duct radius demarcation, and is seen by lower CO values coupled 
with higher C0 2 and NO x values. 

• Jet trajectories may be inferred from charting either the minimal 
temperatures or the maximal 0 2 species concentrations as a function 
of downstream distance from the jet orifice. 

• The temperature and species concentration profiles for the six 
configurations share the same general flow characteristics relative to 
jet penetration dynamics. 

• When jet penetration increases beyond optimal as in the 8-, 9-, 10-, 
and 12-hole cases tested, the jet mass gravitates to and accumulates in 
the central core of the mixer rather than dispersing laterally 
throughout the radius of the mixer. 

• When jet penetration decreases beyond the optimal point as in the 18- 
hole case, a hotter, rich core is allowed to bypass the jet region 
without completing the reaction toward the fuel-lean state. 

• Of the six hole configurations tested, the 14-hole module exhibits jet 
trajectory penetration close to the mixer half-radius and produces the 
most complete reaction and the best uniform mixing by the z/R=2 
plane. 

• The overpenetrating 8-, 9-, and 10-hole cases yielded little reaction 
beyond the z/R=l plane while the optimal 14-hole, the slighty 
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overpenetrating 12-hole, and the slightly underpenetrating 18-hole 
cases showed substantial reaction occurring after the plane. 

• Although the profiles obtained are able to indicate general flow 

* 

characteristics that allow for comparison between mixers, a denser 
data grid would be advantageous in order to capture specific 
gradients more precisely. 

6.2 Recommendations 

The results from this experiment present further questions that may be 
answered by: 

• Acquiring mixture fraction fields using an inert gas tracer to establish 
the extent of mixing versus reaction in the flow field. 

• Obtaining time-resolved measurements of temperature and species 
concentrations to determine the cause of jet flow asymmetry in the 
orifice vicinity. 

• Determining the velocity flow field to quantify mass-flow rates of 
species concentration in the plane. 

• Preheating the mainflow air to investigate the effect of higher inlet 

* 

temperatures on jet reaction and mixing processes as well as on NO x 

9 

production. 

• Comparing mixture and reaction fields from the reacting case to a 
non-reacting case performed at the same momentum-flux ratio J =57. 
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APPENDIX A 


ORIFICE AREA CALCULATION 


Design Conditions : 


mass flow ratio 

MR 

:= 2.5 

main flow cross-sectional area 

^ m 

:= 0.005027- 

jet to main flow momentum-flux ratio 

J := 40 

jet temperature 

T i = 

298K 

main flow temperature 

T m 

:= 1800K 

discharge coefficient 

c d 

= 0.65 


J = (jet density) *(jet velocity^ 

(main density)*(main velocity)2 


where universal gas constant R and pressure P of jets 
and main flow are assumed constant, and velocity is 
represented by: mass flow rate/ density /area. 


Effective Jet Area 


^ jet_eff 


MR- A 


m 


i T i 


J T 


m 


A jet e ff = 808.52* mm" 


Geometric Jet Area 




A: 


jet_geom 


^ jet_eff 


A jet_geom 1-244*10 *mm 
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APPENDIX B 


EMISSION ANALYZER OPERATION PRINCIPLES 


The descriptions of the basic operation principles behind the non- 
dispersive infrared (NDIR), paramagnetic, flame ionization detection (FID), and 
chemiluminescence analyzers were summarized from the Horiba Instruments 
instruction manuals #091891, #091217, and #091216, respectively. 

B.l Infrared Analyzer 

The measurement of CO and C0 2 concentrations in the experiment was 
accomplished by non-dispersive infrared analysis. The non-dispersive infrared 
analyzer (NDIR) (Figure B.l) distinguishes between different molecules by their 
unique infrared absorption band. The absorptivity of the sample in a specific 
band is proportional to the species concentration. 

In the NDIR analyzer, a light source emits infrared light which is 
transformed into intermittent light by a chopper. The intermittent light is then 
passed through a measurement cell where it is absorbed by the sample. The 
absorptivity is compared by the detector to a reading from an adjacent 
reference cell. The change in intensity of light between the sample and 
reference cells causes a membrane in the detector to vibrate, which generates a 
measurable electrical signal. 
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Figure B.l Basic NDIR Detector Components 
(Adapted from Horiba Manual #091891) 


B.2 Flame Ionization Detector (FID) 

Ionization occurs when a hydrocarbon sample is burned with a 
hydrogen flame. In the flame ionization detector (FID) (Figure B.2) a DC 
voltage is applied across two electrodes situated on opposite sides of the flame. 
The electrical potential induces a current, or movement of the ions produced 

# 

from the burned hydrocarbon sample, which is proportional to the number of 
• carbon atoms in the hydrocarbon sample. The current is amplified and 

converted to a voltage differential that is measured. 


ELECTRODES 



Figure B.2 Basic FID Components (Adapted from Horiba Manual #091217) 


B.3 Magneto-Pneumatic Analyzer 

The magneto-pneumatic analyzer (Figure B.3) is used to measure oxygen 
concentrations in a sample by utilizing the paramagnetic properties of gaseous 
oxygen. An uneven magnetic field applied to such a gas causes the molecules 
to migrate toward the portion of the field with the strongest attraction. The 
accumulation of gas raises the pressure at that point. 

The pressure rise AP is directly proportional to the concentration C of 


paramagnetic gas, as shown in equation (B.l): 
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where H represents the magnetic field strength and X represents the magnetic 
susceptibility of the paramagnetic gas. The changes in pressure are converted 
into electrical signals by a condenser microphone. The resulting electrical 
output is linearly proportional to the oxygen concentration. 

To eliminate the pressure rise between sampling, a non-paramagnetic 
gas such as nitrogen (N 2 ) is required in the operation of the analyzer. 


MAGNETIC POLES 


MAGNETIC 
FIELD CELL 



SAMPLE 


CONDENSER 

MICROPHONE 


SAMPLE 


\ V 


AMPLIFIER 


PREAMPLIFIER 


Figure B.3 Basic Paramagnetic Detector Components 
(Adapted from Horiba Manual #091217) 


98 


B.4 Chemiluminescence 

Chemiluminescence is used to detect NO x by initiating a reaction 
sequence between the sample and ozone (0 3 ) to emit light. The sample is first 
routed through a chamber that dissociates the N0 2 present in the sample into 
NO. The gas is then channeled to another chamber where it reacts with 0 3 in 
the following set of reactions listed as equations B.2 and B.3: 


NO + O s -> NO* + 0 2 

(B.2) 

NO/ NO z + ht> 

(B.3) 


where h is Planck's constant, v refers to the frequency of the light emitted, and 
the asterisk designates the electronically excited state of the molecule. The light 
emitted in the second reaction (Equation B.3) is measured by a silicon 
photodiode. The measurement is directly proportional to the concentration of 
the total NO, which includes the N0 2 dissociated prior to reaction with 0 3 , that 
entered the detector reaction chamber. 

Unlike other instruments which measure chemiluminescence with a 
photomultiplier tube, the Horiba model utilizes silicon photodiodes, which 
increases the range of wavelength detection while decreasing noise in the 
reading. 


APPENDIX C 


EXPERIMENTAL CONDITIONS AND SETTINGS 




Atmospheric. 

Room Temperature 


P s'j’D := \ 1.013’ 10 


JN 

2 

m 


T STD = 298K 


Gas Constants 


R C3H8 := °- 1885 ‘ 


kj 

kg-K 


R 


air 


0.2870- 


kj 

kg-K 


Densities of propane 
and air at standard 
conditions . 


PC3H8 = 1-803- kg- m 


-3 


Pa^ = 1.184- kg- m 3 


¥ 


# 


Reacting Crossflow ("main") and Tet Areas and Temperatures 


A main : = 5.027- 10~ 3 -m 2 


T main = 15 00 K 


-^jets_geom 1-237-10 -m 


Tjets = 


Pressure Drop Across Quick-Mixing Orifices 
AP jels : = 448(H>a 

Experimental Mass Flow Rates 

Mass_flow = 0.00296- kg- sec 1 

u n nno7C , _ l Mass flow : ofc = 0.0752- kg- sec 

Mass_flow ^ = 0.0275- kg- sec ~ jets & 

Total rich crossflow mass flow rate (fuel and air): 

Mass_flow ma j ri := Mass_flow a^. + Mass_flow 

Mass_flow = 0.0305- kg- sec 1 
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C.l CALCULATION OF EQUIVALENCE RATIO 


Rich Experimental Fuel-Air Ratio (FAR') : 


FAR 


rich : 


Mass_flow ^ j 

Mass_flow _j r 


FAR rich = 0.108 


i 


¥ 


Stoichiometric Fuel- Air Ratio : 

Stoichiometric reaction assumes no excess fuel or air in products of 
complete combustion. The following is a mass balance of the combustion 
of propane. 

C 3 H8 + 5(O z + 3.71 N 2 ) ~> 3C0 2 + 4H 2 0 + 18.55N 2 


Species Molecular 
Weights 


MW c := 12- 
MW h := 1* 


g 

mol 

g 

mol 


MW 02 := 
MW N2 := 


32- 


g 


mol 


28- 


g 


mol 


MW (23PJ8 3' MW o 8- MW j_j 

MW air := MW 02 3.71- MW N2 


FAR s t 0 i c h 


1- mol- MW C3Hg 
5- mob MW 


FAR s toich 0-065 


* 


* 
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Rich Zone Equivalence Ratio 


^ rich 


FAR r 
FAR stoics 


m 


Overall Fuel- Air Ratio (FAR) : 


FAR overa ii 


Mass_flow 

Mass_flow a i r + Mass_flow j ets 


Overall Equivalence Ratio 


^ overall 


FAR overa u 
FAR stoich 




* 


<t> rich = L66 


FAR overall 0.029 


^ overall ®>45 
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C.2 CALCULATION OF MOMENTUM-FLUX RATIO J 

Calculate main flow momentum flux 


P STD 

P main R <T 

K air 1 main 

4 

P main = 0.235* kg*m 

% 

Mass .flow main 

V - - 

mam o • -A 

e mam mam 

V ma j n = 25.758* m* sec 1 


Mom_flux ma j n P main ^ main M°m_flux m ain 156* kg* m *sec 


Calculate total jet momentum-flux 


P jets = P STD + ^ P jets 

P . 

jets 

p i ets r . t • 

aur 1 jets 

jets 

v . . •= 2 * J 

) ets ' p ■ , 

4 ^ jets 

P jets =0.768* kg* m 3 


Vj ets = 108.022* m* sec _1 

Mom-fluxj ets = P jets' V jets 2 

Mom_fluXj e ^ g =8.96*10 3 *kg*m 1 *sec 2 


Momentum Flux Ratio T 

Mom_flux j ets 
Mom-flux main 

» 

J = 57 


* 


103 


C.3 CALCULATION OF FLOW PARAMETERS 

Reference Velocity 
P STD 

P3 "p . T p o = 1.184*kg*m 3 

K air 1 STD K 3 6 

Mass_flow re f := Mass_flow ma j n + Mass_flow j ets 

Mass_flow re f = 0.106* kg* sec 


^ ref ' 


Mass_flow re £ 
P 3 ^ main 


V re f = 17.742*m*sec 


Mass Flow Ratio MR 

Mass_flow : ets 

MR := — MR = 2.47 

Massjlow main 


Density Ratio DR 
P jets 


DR := 


mam 


Velocity Ratio VR 


V: 


VR 


jets 


V 


DR = 3.26 


mam 


VR = 4.19 
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Effective Tet Area and Discharge Coefficient 


Mass flow 


jets_eff 


jets 


P jets 


AP 


jets 


P jets 


^ jets_eff = 9-062* 10 -4 -m 2 


d 1 


^ jets_eff 
^ jets_geom 


C d = 0.733 




* 


APPENDIX D 
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MASS BALANCE 


Species Molecular Weights 


MW co := 28 

MW 02 := 32 

MW C3H8 ■ 44 

MW C02 := 44 

MW N2 := 28 

MW air := 28.97 

MW NO := 30 

MW c := 12 

MW 44 

MW j^20 := ^ 

MW j_j := 1 



Percent of Element in Air 
Air 02 % := 0.2095 Air jvj2% := 0.7808 


Elemental Stoichiometric Coefficient of Element in H 2 O 
(0^=2 a>Q2 - 0.5 


Experimental Measurement s 


k g 


M fuel : = 0.00296- — 2 - 
sec 


M air * °' 0275 '^ 
sec 


M 


total ^ fuel + ^ air 


M total = 0 - 03, k §* sec 


CO% := 0.13 co 2% ;= °- 0532 

O 2 o /o := 0 HC% •= 0.0129 

NOx% := 0.0000018 

(Assume N2 makes up the rest of the exhaust gas volume.) 

N 2 o /o := 1 - CO% - CO 2 % ~ HC% - NOx% - O 2 o /o N 2 o /o 


0.804 
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Calculate Mass Flow Rate of each element entering system: 


M H : = M fuel' MW H- w 

MW C3H8 

= 5.382*10 4 • kg* sec 

3 

-1 

M C =M M -MW f 

MW C3H8 

M q = 0.002* kg* sec 

Air 02% 

, -i 

M 02 =M alr MW 02 - 

iVAVV air 

Mq 2 = 0.006* kg* sec 

Air N2% 

, -l 

M N2 =M air MW N2 - 

iVAVV air 

Mn 2 = 0.021*kg*sec 


Calculate Mass Flow of Dry Exhaust Gas Exiting System: 

MW Gd = (MW co CO%) + (MW C02 C0 2 %) + (MW N q-NOx%) + (MW H( ~'HC%) + (MVV q 2 0 + (MW n 2 ' n 2%) 

MW Gd = 29 - 058 
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Calculate Mole Fraction of Water per Element and Compare: 


Carbon Balance 

MW C -M tota] - (c«tCO 2 % tHC% ) - M C MW Gd 
VmOC MW C M total - (CO% -hCO 2% + HC%) - M c (MW Gd - MW H20 ) 


ymoc 1 


Hydrogen Balance 

mw h m tota] - (2.67hc%) - m h mw 


y H20H 


Gd 


MW H'M total' (2.67-HC%- (0 J - M H - (mw Gd - MW H2Q ) 


y H20H : 


Oxygen Balance 


y H200 ' 


MW 02 M total - (0.5-CO%-hCO 2 o /o + 0.5 nOx%+O 2% ) - m Q2 mw 


Gd 


MW 02 M total' (0.5-CO%d-CO 2 o /o -h0.5-NOx%-hO 2 o /o - <0 02) ~ M 02" Gd ” MW H2o) 


y H200 : 


Nitrogen Balance 


y H20N 


MW N 2 -M total- (0-5-NQx%+N 2 . /o ) - M N2 -MW Gd 


MW M2‘M total (0.5-NOx%-|-N 2 „ /o ) - M N2 - (MW Gd - MW h2 q) 


y H20N = 




0.029 


0.222 


0.158 


0.181 


APPENDIX E 


TEMPERATURE AND SPECIES CONCENTRATION HISTOGRAMS 


Histogram plots of temperature and species concentration for the 8-hole 
module were shown previously in Figures 5.1, 5.2a, and 5.2b. For comparative 
purposes, the raw data histograms for the other module cases (9-, 10-, 12-, 14-, 
and 18-hole) are presented on the following pages. 

Figures E.la, E.lb show the temperature histograms for these cases. 
Each column of histograms corresponds to the particular module case, while 
each row represents the plane at which the 16 points of data were obtained. 

The species concentration profiles are shown in Figures E.2a, E.2b for the 
9-hole case; Figures E.3a, E.3b for the 10-hole case; Figures E.4a, E.4b for the 12- 
hole case; Figures E.5a, E.5b for the 14-hole case; and Figures E.6a, E.6b for the 
18-hole case. Each column corresponds to the measured specie (either 0 2 , C0 2 , 
CO, HC, or NO x ), and each row represents the plane that the data were 
obtained. 
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Figure E.3a 0 2 , C0 2 , CO Concentration Distribution for the 10-Hole Module 
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Figure E.5a 0 2 , C0 2 , CO Concentration Distribution for the 14-Hole Module 
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Figure E.5b HC and NO x Concentration Distribution for the 14-Hole Module 
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Figure E.6a 0 2 , C0 2> CO Concentration Distribution for the 18-Hole Module 














120 


Plane 6 


QzOzQ: z 


Plane 1 


ttt 


HC(%) 


2 

1.5 

1 

0.5 

0 


rflinsoirnio 


2 

1.5 

1 

0.5 

0 


. 


- 


T-comrv-OT-com 


2 

1.5 

1 

0.5 


m 


0 f i* il —4 — JU- -- f— J t -it -j 


T-coior^oiT- com 


2 

1.5 

1 

0.5 

0 


T-nwso5i-oin 


NO x (ppm) 



36 j 
27 - 

■ lull 


18 - 

mrnaniii 


9 - 

iii ii mini in 


0-1 

ill mi i I i i ill 1 1 1 


t- co in is. 


36 

27 

18 

9 

0 


, 1 fi | | 


IK II 

miifffMiif 


-»-coLor^05T-com 


36 

27 

18 

9 

0 


EMJUL&JL 


lihll.i.1 


T-cowh-o)T-com 



36 1 
27 - 

: ’ . 

n _ . . ... . 

18 - 

' 

i I I i m i ram 

9 . 


Mil 11 11 11 ill 

0 - 

l:ljB t aJ.fl.a JLeuII 


i-coioi^cD-r-com 


Plane 6 
z/R=2 


Plane 5 
z/R=l 


Plane 4 
z/R=d/R 


Plane 3 
z/R=(d/2)/R 


4 

A 


R 



2 

1.5 

1 

0.5 

0 


T-cowh«.<x>T-com 


2 

1.5 

1 

0.5 

0 



36 -j 
27 

18 

Q 


ppgMliiigaiilil 


Tijjijiiiiprp 


i iililiiR imiK 


S H JIjajJl.fl.flfl B gallon 


i-oiflso>r-com 



36 
27 
18 
9 . 


(IIIII 

a a b fi 1 

III 

III 

II 

’Vi;' • • .1 

111 111111 Mil III! 

0 

joLB.ft.H j.b a a li oi.fli.B.a 


Plane 2 
z/R=0 


Plane 1 
z/R= -1 


4 


POINT 


POINT 


Figure E.6b HC and NO x Concentration Distribution for the 18-Hole Module 
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